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We consider the generation of primordial curvature perturbation by general non-Abelian vector 
fields without committing to a particular group. Self-interactions of non-Abelian fields make the 
field perturbation non-Gaussian. We calculate the bispectrum of the field perturbation using the in- 
in formalism at tree level. The bispectrum is dominated by the classical evolution of fields outside 
- - - the horizon. In view of this we show that the dominant contribution can be obtained from the 

(N ■ homogeneous classical equation of motion. Then we calculate the power spectrum of the curvature 

. perturbation. The anisotropy in spectrum is suppressed by the number of fields. This makes it 

' possible for vector fields to be responsible for the total curvature perturbation in the Universe 

without violating observational bounds on statistical anisotropy. The bispectrum of the curvature 
(~| , perturbation is also anisotropic. Finally we give an example of the end-of-inflation scenario in which 

' the curvature perturbation is generated by vector gauge fields through varying gauge coupling 

' constant(s), which in covariant derivatives couples the Higgs field to the vector fields. We find that 

, reasonably large gauge groups may result in the observable anisotropy in the power spectrum of the 

curvature perturbation. 

O: 
^: 

^ • I. INTRODUCTION 

Q ' Inflation was proposed to alleviate the horizon and flatness problems of the Hot Big Bang cosmology [JQl- Currently 
^ it still is arguably the most compelling mechanism to explain the high degree of homogeneity and isotropy of the 
c/3 Universe as well as its flatness. The case for inflation was further strengthened after the release of first detailed 
observations of the Cosmic Microwave Background (CMB) radiation. The spectrum of the temperature perturbation 
of the CMB was found to be consistent with infiationary predictions and ruled out the rival theory of cosmic strings as 
, the primary origin of this perturbation. Moreover, with the increasing precision of observational data it is becoming 
^ possible not only to falsify the competing theories of the origin of the primordial perturbation, but to falsify different 
models of inflation themselves. Measurements of temperature irregularities in the CMB sky provide a powerful tool 
to probe the physics of the very early Universe and the increasing precision allows to do it in a more and more detail. 
[ The simplest inflationary models, with a scalar field driving infiation and producing the primordial curvature 
perturbation C, are still consistent with current CMB data. But several tentative anomalies, which are persistently 
found in all WMAP data releases, might suggest a need of more complex models. For example, the observed power 
asymmetry the alignment of low-/ CMB multipoles or a deep cold spot in the southern Galactic hemisphere 

[9, 10]. If the origin of these anomalies is confirmed to be primordial it will imply some degree of the statistical 
inhomogeneity and/or anisotropy of the primordial curvature perturbation. Such anomalies cannot be explained by 
the simplest infiationary models invoking only scalar fields. Scalar fields do not generate statistical anisotropy as they 
do not choose a preferred direction. But vector fields do, and if statistical anisotropy is established to be of primordial 
$^ origin, the most natural way to explain it is with effects of vector fields. 

There are two ways through which non-negligible contribution of vector fields to the evolution of the universe can 
induce statistical anisotropy. First, if a vector field have an effect on the global expansion of the universe, the latter 
will be anisotropic. The anisotropic expansion during inflation causes statistically anisotropic quantum fluctuations 
of light scalar fields. When these fiuctuations become classical and cause perturbations in the metric, the latter are 
statistically anisotropic too as well as the temperature irregularities of CMB. This mechanism was first considered 
in fll'|. On the other hand, the energy density of vector fields can be negligible, so that the global expansion is 
approximately isotropic, but perturbations of vector fields themselves generate or contribute to (.■ In general quantum 
fluctuations of vector fields are statistically anisotropic making their contribution to C statistically anisotropic too. 
The generation of C by vector fields was first considered in Ref. [l2|. However, this and several subsequent papers 
[T3I . [iJI did not consider statistical anisotropy, which was first considered in Ref. [l^ . The comprehensive study of 
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statistical anisotropy generated by a vector field can be found in Ref. fld\ (see also Ref. [T3| for computations of 
anisotropic non-linearity parameter /nl)- 

Recently the interest in statistical anisotropy and vector fields has grown considerably. Authors of Refs. (l8l - [23| 
studied two-, three- and four-point correlators of the curvature perturbation generated by vector fields in detail. While 
in Refs. [2^ . [25| a model with a massive vector curvaton field is presented in which the vector field can generate both 
statistically isotropic and anisotropic perturbation. In these works a negligible contribution of vector fields to the 
global expansion rate is assumed. 

Another line of research was concentrated on effects of anisotropic inflation on the curvature perturbation (26l-l28| . 
The universe during inflation expands anisotropically if the backreaction of the vector field is non-negligible. Such 
setup is considered in Refs. [Til |29| - |3^ . Particularly interesting are the results of Refs. [13, [13 ■ Authors of these 
papers considered a model with time varying kinetic function of an Abelian vector field of the form f {t) F^^F^^'^ , 
where is the field strength tensor F^i, = d^A,j — d,jA^. It was shown that if / {t) is modulated by the inflaton, 
then the scaling of the form / (x a"'* (a being a scale factor) is an attractor solution for a large parameter space. This 
is very significant as such scaling leads to the flat perturbation spectrum for the vector field [3, [Hi • fF^ coupling 
also induces anisotropy in inflationary expansion of the order of the slow- roll parameter. In addition in Ref. [ss"] it was 
shown that the vector field backreaction slows down the inflaton, i.e. the inflaton potential is effectively "flattened". 
The backreaction of non- Abelian vector fields with the time- varying kinetic function was also considered in Ref. [1^ . 

The interest in detecting statistical anisotropy in the CMB is increasing too [36l - [4l] |. The anisotropy in the power 
spectrum of the primordial curvature perturbation can be parametrized as 



(fc) = V'^° (fc) 



1 + .gc (fc) ( n • k ■ 



(1) 



where only a quadrupole term is kept, n and k are unit vectors and the amplitude (k) is in general a function 
of wavenumber k. First results of measuring were given in Ref. (37| . After correcting a mistake in this work 
Refs. (38l.[39| published consistent results detecting the departure from statistical isotropy at very high 9cr significance 
level = 0.29 ±0.031. They also found that the preferred direction n is very close to the ecliptic pole {I, h) = (96, 30). 
The proximity of n to the ecliptic pole suggests very strongly that the detected is due to a systematic effect. This was 
indeed discussed in Refs. (39l . [40| were they investigated whether a non-zero g^ could be accounted for by the WMAP 
beam asymmetry or other systematics. Unfortunately, both works reached somewhat contradicting conclusions. The 
authors of Ref. [S^l could not determine a systematic effect causing such a large , while Ref. [l^l claims it is due to the 
beam asymmetry and give the bound \g(^\ < 0.07. Despite this, as the precision of measurements increase, the prospect 
of detecting statistical anisotropy of primordial origin is very exciting. Such anisotropy offers a new observable to 
probe the physics of the very early Universe. Indeed, as shown in Refs. (sgI . |41] | with the Planck temperature data 
alone it will be possible to constrain g(^ with an accuracy of 0.01 (2a). In addition, the polarization data alone will 
offer 0.03 accuracy as the consistency check. Furthermore, an extended Planck mission can constrain the spectral 
index of g^ oc fc"^ to an accuracy of Aq ^ 0.3 (Icr). 



The interest in measuring the anisotropy of higher order correlators is increasin g to o 



This is justified. 



because higher order correlators can be predominantly anisotropic as shown in Refs. [I^-[l9, 21, 22, [^[l^. This was 
specifically emphasized in Refs. (TtI. [23. [25l [45| were non-linearity parameter /nl was calculated for several models. It 
was found that even if anisotropy in the power spectrum is subdominant, /nl can be predominantly anisotropic with 
the same preferred direction fi as the spectrum. Moreover, the magnitude of /nl is proportional to anisotropy in the 
spectrum g^. If such non-Gaussianity is detected it would be a smoking gun for the vector field contribution to the 
primordial curvature per turbation. Impacts of statistically anisotropic primordial perturbation on CMB observables 
were studied in Refs. |28l [46|. 

Instead of vector fields being responsible only for the curvature perturbation, Refs. (ZtI - Is^ also propose models in 
which vector fields drive inflation. Such models face a challenge of making expansion of the universe predominantly 
isotropic. The two suggestions proposed are either introducing an orthogonal triad of vector fields or a large number 
of randomly oriented ones. Particularly, authors of Refs. ^Sll . [s^ also consider non- Abelian vector fields. Their setup 
consists of an orthogonal triad of equal norm vector fields of SU (2) group. In this case the total energy-momentum 
tensor of vector fields is isotropic and if they dominate the universe, the latter infiates isotropically. 

The study of vector fields in inflationary cosmology are not only interesting from the phenomenological point of 
view, as means of explaining CMB anomalies. Such models are also interesting from the theoretical perspective. The 
possibility of vector flelds affecting or generating the total curvature perturbation in the Universe opens a new window 
for inflationary model building. We might no longer need direct involvement of scalars to create it could be created 
by vector fields with, for example, varying couplings. In particle physics models non- Abelian vector fields are much 
more common than Abelian ones. In addition, as large non-Abelian gauge groups have many vector fields, assuming 
their random orientation, it is natural to expect the suppression of statistical anisotropy of (. 
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In this paper we study the generation of C by non-Abehan vector fields. The effects of SU (2) non-Abehan vector 
fields on the curvature perturbation were first studied in Refs. ^2Tl-|23j. In this paper we demonstrate that non-Gaussian 
correlators are dominated by the interactions outside horizon. First, we make use of the full quantum formalism, the 
so called "in-in formalism", to calculate the bispectrum at tree level and show that it is dominated by the classical part. 
That is, by interactions of fields after horizon crossing. Since this is the case, correlation functions can be calculated 
much easier using only the homogeneous equation of motion of vector fields. We perform this calculation explicitly 
and show that the result obtained is the same as the dominant part of the result calculated in the in-in formalism. 

We also calculate the spectrum and bispectrum of the curvature perturbation. It is found that the anisotropy in 
the spectrum of C, is suppressed by the number of vector fields involved in generating C, (assuming these are oriented 
randomly). Thus with a large enough gauge group C can be generated solely from vector fields without violating 
observational bound on the statistical anisotropy. However, it is possible that a small detectable anisotropy remains. 

The bispectrum of C, from non-Abelian vector fields is also anisotropic. Although the form of anisotropy is more 
complicated than in the single field case as it involves not one but many preferred directions. 

In the last section we give a simple example of a non-Abelian gauge field generating C,. We use the end-of-inflation 
scenario in which C, is generated by the gauge fields through varying gauge coupling constant (s) in the covariant 
derivative, which couples Higgs field to the gauge fields. 

The paper is organized as follows. In section|ll]we define the Lagrangian and the setup of our model. In section ITTTl 
the bispectrum at tree level is calculated using the quantum in-in formalism. In section IIVI it is shown that the 
classical calculation of the bispectrum from the homogeneous equation of motion gives exactly the same result as the 
dominant part of the bispectrum in the previous section. Using the 5N formula, the correlators of C, are calculated 
in section [V] Finally, an example of a mechanism for the generation of C from non-Abelian gauge fields is given in 
section IVTl It is generated through the gauge couplings of the Higgs field. A summary is presented in section rVIII 

In this paper we use natural units, where c = h = \ and Newton's gravitational constant is SttG = rn^^ ^ where mpi 
is the reduced Planck mass. 



II. THE LAGRANGIAN 
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Consider a general Lagrangian of non-Abelian vector fields 

1 

where the field strength tensor is 

Fl,=d^Al~d,Al + i'^^^A\A'l. (3) 

We do not specify the gauge group, f''^ are structure constants of the Lie algebra of any non-Abelian group and they 
are antisymmetric in permutations of indices a, h and c. 

In this paper we are interested in a time varying g but we do not specify the origin of this variation. The modulation 
of g might be due to some scalar degree of freedom. The infiaton itself can modulate g (for such models with Abelian 
vector field see Refs. [29l l3ll - [33| and non-Abelian ones see Ref. [35l|), but in this paper we do not need to specify the 
origin of time dependence of g. 

Let us recast the Lagrangian in Eq. ([2]) by factoring g (£) into a constant and time dependent parts in the following 
way 



5(0^5c/VM- (4) 

(/c is the constant value of g when the modulating degree of freedom is stabilized at time t^. At this moment the 

kinetic function f (t=) = 1. The constant can be absorbed into the field strength tensor by the field redefinition 
Ji^ = A°/5c- Dropping tildes the Lagrangian in Eq. ^ then becomes 

£^-i/F;,i^r, (5) 

where the redefined field strength tensor becomes 

f;, = d^Al - d^Al + ffcf^'^^A^ (6) 



with g^ being a self-coupling constant. 
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Using the gauge freedom of the non-Abehan vector field in Eq. ^ we can always choose a gauge in which any 
given space-time component of all the vector fields is zero [s^- However, in general it is not possible to make all four 
components of all vector fields in a group vanish by a gauge choice. This is in contrast to the Abelian field, where one 
can always choose a gauge in which the classical or homogeneous part of the field is zero. For the rest of the paper let 
us choose the vanishing component to be the temporal one so that Aq = for all a. For the remainder of the paper 
we use the convention of space indices denoted by subscripts and gauge ones by superscripts except in those cases 
where equations are written in Lorentz covariant four-vector form. 

The energy-momentum tensor of the massless vector field is anisotropic [l^ . If the energy density of such a field is 
non-negligible the expansion of the universe becomes anisotropic. To avoid excessive large scale anisotropy one can 
introduce a large number of vector fields, suppressing the aniso trop y by a factor of Vf^, where Af is the number of 
vector fields. Such mechanism is employed in vector inflation [48|. In our setup this could be achieved by taking 
a very large gauge group. Another possibility is to introduce three identical, orthogonal vector fields [4^ [SJI- This 
option was recently explored with SU (2) vector fields in Refs. [sills^l- In this paper we consider a third possibility, 
namely we assume that the energy density of the vector field is negligible during inflation. In other words, we neglect a 
backreaction of the vector field on the expansion of the universe. However vector fields can still generate the curvature 
perturbation. This might happen, for example, in the vector curvaton [12] or end-of-inflation [l5| scenarios. In the 
former case, the non- Abelian vector fields must acquire a mass through a Higgs mechanism prior to generating C. An 
example of the latter case will be given in section IVII 

Taking the contribution of the vector field to the total energy density to be negligible, inflation can be assumed 
to be isotropic and we can use the Friedmann-Lemaitre-Robertson- Walker (FLRW) background with the metric 
5/11/ = diag [l, —a^ (t) , —a^ (t) , — (t)] , where t is the cosmic time. We will also use the conformal time t = J dt/a (t) 
in the paper. 

Expanding the Lagrangian in Eq. ([5]) in the FLRW background we flnd a term 

in which the scale factor appears explicitly. But the normalization of the scale factor is arbitrary while the Lagrangian 
is a physical quantity related to the energy of the system and cannot contain arbitrary normalizable factors. The 
appearance of a in Eq. ([7]) is due to the fact that the vector fleld A"; is defined with respect to the comoving coordinates. 
The physical vector field, defined with respect to the physical coordinates is Af/a [l2l [l^ [isj. Therefore, it will be 
useful to define a physical, canonically normalized vector field 

= Vf^. (8) 

a 

We also use Fourier modes JVKf (k) of the perturbation of W"" 

/aSi. 
—-^5Wt{i^)e^''-. (9) 
(ztt) 



III. CORRELATION FUNCTIONS IN THE IN-IN FORMALISM 



A. The in-in Formalism 

In this section we calculate the three point correlation function of vector field perturbation. Following Refs. (55l.[56| 
we use the so called "in-in formalism". In our setup, the calculation is simplified due to the assumption that the energy 
density of vector fields during inflation is negligible, which means that it is not necessary to consider the corresponding 
perturbation of the metric. 

In the in-in formalism expectation values (0| O |0) are calculated in the interaction picture. In this picture both, 
fleld operators and state vectors, varies with time. The variation of the former is governed by the equation of motion 
of the free field, while the latter evolves due to interaction terms. To find free field operators for non- Abelian vector 
fields we note that each massless vector field has two degrees of freedom. The field is quantized by promoting these 
degrees of freedom to operators with canonical commutation relations. Thus SW°' (k) in Eq. ^ becomes 

SW,^k,r)= J2 [ef (k)«^(fc,r)aX(k)-ef (-k)z«*(fc,r)af (-k)] , (10) 

A=L,R 
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where r is the conformal time. The raising and lowering operators in this equation satisfy canonical commutation 
relations 



al (k) , a\\ (-k') = (27r)^ S (k + k') S^yS^b 



(11) 



with others being zero, while the left- and right-handed circular polarization vectors e^ are chosen in such a way that 



(-) 



k ) . With k = (0, 0, k) they become 

ei^(k) =i=(l,z,0) and ef 



(k)=i=(l,-.0). 



(12) 



In the interaction picture SW^ (k, r) is a free quantum field. Thus the conjugate pair {w, w*} in Eq. pop are solutions 
of a free vector field equation of motion. In principle {w,w*} should have polarization and group indices, e.g. w^. 
However, the Lagrangian in Eq. ^ does not have parity violating terms and both polarization modes satisfy the 
same equation of motion. The same is true for the gauge index, all vector fields satisfy the same free field equation 
of motion with the same initial conditions. Thus to simplify the notation we dropped polarization and gauge indices 
out. 

In Refs. [U, [2^ it was found that w satisfies the same equation of motion as the scalar field if the kinetic function 
during inflation varies with time as / cx a^^^^. For a massless field in de Sitter background this equation becomes 

2 



zij -I- 3Hw 



-I w = 0, 

a. 



(13) 



where dots denote derivatives with respect to the cosmic time t. The solution of this equation is very well known. 
With Bunch-Davies vacuum initial conditions it is given in conformal time by 



H 

TIP 



(1 — ikr) e 



With these definitions the Wightman function becomes 



(2^)^ 6 (k + k') SabT^ (k) w (fe, r) w* (fc, r') , 



(14) 



(15) 



where the tensor ^k^ is defined by (igI . l45j 

T^(k)-ef(k)ef(k)+ef (k)e,^(k)=^ 

N-point correlation functions of the vector field perturbation are calculated in the interaction picture as vacuum 
expectation values of the form 



(16) 



gN (Xl,X2, . . . ,XAr) 



U-^SW/ (xi, r) SW^ (X2, r) . . . SWj^ (xat, t) U 



, 



(17) 



where, to simplify notation, we also suppressed gauge and space indices in the function g^. The unitary operator U 
is given by 

f7 = exp|-z^ Hi„t(r')dr'| (18) 

with To being some early time when the mode of interest is deep within the horizon. The interaction Hamiltonian 
Hint can be found from the interaction terms of the Lagrangian in Eq. ([5]). With the temporal gauge (v4g — 0) these 
terms become 



1 



gcf'" {diA"^ ~ djAf) A\A] + -gli'''"'i'"^''A\A]AiA] 



(19) 



In this paper we calculate the three point correlation function at tree level, thus only the third order interaction 
Hamiltonian is considered. For the physical, canonically normalized field it is i?int = ii'^t ^ ^'^'^^ 



H\J^, where 



U'>^^ = a-" (r) / d^x ^i'''"'dm''&Wtm. 



■9c cabc ' 



(4) 



V7 

a" (r) / d^x ^pfccpde ^ pdcp&e^ wlm'^mfm. 



(20) 
(21) 
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The factor in these equations is due to -^Z— det [g^i/]. 



As it is clear from Eqs. ([20ll and (|2T|) g (t) = gc/\/fW) Eq. (g]) is the strength of self-coupling for the canonically 
normalized vector field. To keep quantum calculations under control this coupling must be ensured to be small. It must 
also be small for quantum fluctuations of interacting fields to become classical after horizon exit |57|. Furthermore, the 
perturbation spectrum of massless non-interacting vector fields is flat if the kinetic function is of the form / cx a~^^^ 
pi. [25|. In order to preserve approximate flatness of the perturbation spectrum, interaction terms must be small. 
Because the variation of / (i) is exponential at a time t < ts with / (tg) = 1 by deflnition, g ~ gc/V7 is small only 
when / c>c a~^. We will assume this to be the case for the rest of the paper. 

With the weak self-coupling Eq. can be expanded in powers of gc/\/J- As we are interested in the tree level 
contribution to the bispectrum it is enough to keep only the flrst order term in Eq. (|17p after which it becomes (sgI [58| 



53 



dr'/o 



SW/ (xi, r) 6Wl (X2, t) 6W!: (xg, r) , H,^^ (r') 



(22) 



In what follows g^'' will denote the three point function with the Hamiltonian in Eq. (|20p and 53'*'' in Eq. (PTjl . so 



that the total function is g^ 



B. The Three Point Correlation Function from the Quartic Term 



Let us consider flrst the ^g*'' term. Going to the momentum space we find 



g^''^ (ki,k2,k3) = (^f'^bc^'^de ^ ^adc^abe^ 



rb9c 

fo 



2 /-Tond 



dr'flS (r') 



{2nf 



S (qi + q2 + qs) X 



-i(0 



SW/ (ki, r) SWi (k2, r) SW^ (kg, r) SWf (qi, r') SWf (q2, r') SWf (qg, r') o) (^3) 



xRc 

where Re [. . .] denotes the real part and Tend is the conformal time at the end of inflation. In this expression we used 
the fact that Wj' is slowly rolling due to the smallness of self coupling term (more on this in section IIVB[) . We have 
also used / — /oa~^, where /o is some initial value. The correlator can be evaluated using Wick's theorem. With 
Eq. (fT5|) after some tedious algebra we flnd 



gi'^ (ki, k2, kg) = - {27rf S (ki + k2 + kg) ,^—7<^^^^^ (kl, ks, ks) (fcl, fc2, k 



(24) 



The function Ti^^^'^ depends only on the direction of three vectors (ki,k2,kg) and thus quantifles the anisotropy of 



the three point correlator function. The full expression is given by 

7;L'i^^" (kl,k2,k3) ^ WtT^^ (ki) r,f^. (kg) {i-bh^-9f ^^agH^a. 



+W('T^j (k2) T^j (ks) (f^^f p/'» + p/gpb/.^ 

+w^te (ki) T:,fy (k2) {r'H'^''^ + f^'^siabf^ 



(25) 



where (k) is defined in Eq. (|16p . I^^'^ is the integral of the form 



= ^Re 
Jo 



dr' (1 - ikir') (1 - ik2T') (1 - i/cgr') e**^'^' 



(26) 



with fct = fcl + fc2 + fcs- This integral is calculated in Appendix [Xj Assuming all three fc's crosses the horizon at a 
similar time it is equal to 



r(4) 



h 4! 



6e*^'= 



Ki+K2]+ 26^^" Ki - iK; 



300 I ^ ^ 49 y 



(27) 
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where 

and K.^l^. (28) 

7 « 0.577 is Euler-Mascheroni constant and Nk = — In (|fcfTcnd|) is the number of e-folds from when kt exits the horizon 
to the end of inflation. For the cosmological scales N/^ ~ 60. We recognize the dominant term to be proportional 
to the dominant term of the three point correlation functions in Refs. pSi. [59|. This term is the contribution to the 
correlation function from the superhorizon evolution of the fields. The set up in our case is somewhat different. The 
authors of Refs. [58, 59] considered a field with constant strength of self-coupling, while in our case the self-coupling of 
canonically normalized fields is varying with time, g — gc/ \fj oc c? . From Eq. ([27| we see that this variation enhances 
additional modes, both, when they are created at the horizon exit and during the evolution of the field outside the 
horizon. However, these modes are subdominant. Taking only the dominant contribution to g^^ we find 

gf = - (2.)^ 8 (k, + k, + k3) r^^"^ (kr, k„ ka) ^ (29) 

To evaluate this equation we used /end ~ fk exp (— 4A^fc), where fk and /end are the values of / at the horizon crossing 
and the end of inflation respectively. 



C. The Three Point Correlation Function from the Cubic Term 



From Eq. ([20|) we expect that the three point correlation function from the cubic term g^^^ is suppressed by a factor 
of p compared to g^^\ where p — k/a is the modulus of the physical momentum. In this section we show that this is 
indeed the case. 173'^^ is calculated along the same lines as g^^\ Taking a Fourier transform of Eq. ([20| from Eq. ((22|) 



we find 



(3) 
53 



9c _pfcc 



Vfo 
X 2i Im 



dr' (r') 



S (qi + q2 + qs) qu x 



SW/ (ki, r) SWl (k2, r) SW^^ (kg, r) SW^ (qi, r') SWj^ (qa, r') SW^ (qg, r') 
where Im [. . .] denotes the imaginary part. Using Wick's theorem and Eq. ([T5|) we calculate 

gi'^ = - (2^)^ S (ki +k,+ k3) -^7^1^^'^ (ki, k2, k3) (fci, ks) , 

where the anisotropy of the three point correlation function is given by 



7;L'i/^nki,k„k3 



T^y (ki) T^, (ka) r/„ (kg) (pi - pa), + 
+Tl^ (ki) Tf„, (ka) T^i (ka) (P2 - Pi), + 



+TS (ki) T^^ (k2) 7f„ (kg) (P3 - P2), 
p = k/flcnd in this expression is the physical momentum, Oond is the scale factor at the end of inflation and 



/(3) ^ 



=Im 



dr' (r') (1 - ikir') (1 - ifcar') (I - ikjr') e^'^''^' 



Using the method explained in Appendix [A] it is calculated to be 



/(3) 



--K1 + K2 



l[\-K,+2K, 



(30) 



(31) 



(32) 



(33) 



(34) 



The first, dominant term, is due to the evolution after horizon exit. Neglecting the subdominant term becomes 



(3) 
53 



I {2nf 6 (ki + k2 + k3) Ti^±^'' (ki, k2, k3) 



(35) 
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The anisotropic term T'^'*-' from quartic interactions in Eq. ([55]) is proportional to the homogeneous part of the vector 
field W > H,^ while T^^'> from cubic interactions is proportional to the physical momentum p, which for cosmological 

scales are p H. Thus gi^'' ^ g^^ and the dominant contribution to the three point correlation function is from 



quartic terms. Moreover, the dominant contribution in g'^^ itself is from the classical evolution of fields. 

IV. BISPECTRUM FROM THE CLASSICAL EVOLUTION 

A. Classicality 

In the last section we used perturbative quantum field theory to calculate correlators of the field perturbation. As 
the results in Eqs. (|29p and ([55t show those correlators are dominated by the interaction of fields after horizon exit. 
Furthermore, as correlators with derivative couplings are suppressed by the factor fc/a <C H, it suggests that we can 
obtain the dominant contribution to correlators by a simpler method: from the classical equation of motion of the 
homogeneous field. 

It is well known that after a mode k of a free light quantum field crosses the horizon, i.e. when k/aH — !■ 0, the 
phase of the mode function becomes constant and field operator in the Heisenberg picture can be written in a form 

Xk (t) ^Sk + SLk) I where Xk is made real by an arbitrary phase rotation [60l [6l|. In this limit all commutators of 
fields vanish and the eigenvector of a field operator at some particular time remains an eigenvector thereafter. This is 
a cosmological analogue of quantum decoherence. When this happens, quantum fields are well described by classical 
stochastic functions and we say that the field enters into 'classical evolution'. This is true for light free quantum fields. 
Canonical massless vector fields, however, do not become classical as their Lagrangian is invariant under conformal 
transformation to flat space-time. This is a case, for example, with U {!) vector fleld with minimal kinetic term [(3^. 
In our case, conformal invariance of the vector field is broken by the time varying kinetic function / in Eq. But 
the ques tion remains whether self-interaction terms do not prevent the non-Abelian field from becoming classical. In 
Ref. 571 it was shown that after horizon crossing the interacting field does become classical if the interaction is weak; 
specifically if U in Eq. (|18p is sufficiently close to unity, C/ ~ 1. As was discussed before this can only happen if / is 
a decreasing function in time. With a constraint of the fiat perturbation spectrum this means / cx 

In this section we show that correlator functions of the non-Abelian vector field perturbation can be calculated 
using the classical equation of motion. 

B. The Equation of Motion and the Power Spectrum 

Extremising the action with the Lagrangian in Eq. ([5]) we obtain the field equation for non-Abelian vector fields 

5A + aAln y/-det[5^.] (/^^^) - /5cf"''"4^^r = 0, (36) 
Taking the spatial component of this equation k = i and adopting a temporal gauge Aq — we get 

A'^ +(^H+J^ - {9jdjA'^ - d.djA'j) - a-^gj""^'' [2 {A^jAf + d.A^A^^) - gJ'''^''A''^A^A'r] = 0. (37) 

We are interested in superhorizon evolution of the vector field perturbation. On these scales derivative terms are 
negligible and Eq. (|37p can be written as 

A'l + Ih + j:j A'l + glr''^i'"^''a-'^A]AjAl = 0, (38) 



^ For the perturbative approach to be valid W > &W must hold. The typical value of the field perturbation is 5W ^ ~ H/2it, 

resulting \nW>H. 
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which is the same as that of the homogeneous mode of the field. Changing to the physical canonically normalized 
vector field in Eq. ([5]), with / oc this equation is transformed into 

Wj" + SHWl" + !^i''''''i'"^'=W^WfW,^ = 0, (39) 
/ 

which is reminiscent of an interacting scalar field. As fields in Eq. (I39p are intended to generate the curvature 
perturbation, they must retain an approximate scale invariance in accordance with observations. For this to be the 
case W.l^ must be almost a free field, in other words W/* must be rolling slowly. In analogy to Ref. [6^ we introduce 
slow-roll parameters 



^ ' and Vt^:^, (40) 



~ 2mpiiJ ~ 3i?2 



where = dV^'/dW^ and 



.,2 



J 

and require l^fjl < 1 and \rifj'\ ~|j(^)^<f. As structure constants are of order unity, the slow- roll conditions 

mean that the strength of self-coupling g {t) = ffc/v/W Eq. ([4]) of canonically normalized field is small, i.e. 
9c/ f 1. This is easily achieved when cosmological scales exit the horizon. Because / cx a^* is an exponentially 
decaying function with / (tg) = 1, the self-coupling g^/f is exponentially suppressed. Note however that although 
^ 1 is easily satisfied when cosmological scales exit the horizon, this condition must hold up until tg. Even if 
the evolution becomes strongly non-linear after cosmological scales crosses the horizon, all scales are affected. Thus 
not only g^// must be small at horizon crossing but it must remain small when / {ts) = 1, i.e. g^ < 1. We assume 
this to be the case. 

Following Ref. [s^ we decompose the vector field as 

= + SW^ + ^6W^ + ..., (42) 

where the field notation without space indices means the modulus, e.g. = The first term in this expression 

is the homogeneous field. For the rest of the paper we will have no use of the total vector field W"", thus we drop 
out the subscript '0' from the homogeneous mode and denote it simply by W"". The second term in Eq. (|42|) is the 
perturbation and later terms are higher orders in (5W^". This expansion is not unique and to determine 6Wi some 
auxiliary conditions need to be imposed [s^ HJ- We choose 6Wi in such a way that its equation of motion is linear, 
i.e. Eq. (|39p without the last term. It follows then that 5Wi obeys the Gaussian statistics and its two point correlation 
function is 

(<5M/f, (k) SW',^ (k')) = {2n f S (k + k') ^d^tT^^ (k) (k) , (43) 

where ^k^ is defined in Eq. and the Fourier transform of SW^^ is defined in Eq. ^ . We also used the fact that 

there is no correlation between left- and right-handed modes, hence the Kronecker delta Sat- In de Sitter inflation 
the power spectrum is scale invariant = (7J/27r)^. In the slow-roll inflation the spectrum acquires weak scale 
dependence due to slowly increasing horizon size. Assuming approximately constant H / it is 

hV { k 



The spectral tilt '— 2e' is due to the slight increase of the Hubble horizon during slow-roll inflation, which is 
parametrized by e = —H/H^. Because e > 0, each subsequent k mode crosses a horizon of larger size making 
the amplitude of perturbation smaller. 

The spectrum in Eq. (|44p is for the non-interacting part 6Wi of the field perturbation. However, the total power 
spectrum will have additional scale dependence 
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where \ri\ = lilfj] is the modulus of the slow-roh parameter matrix in Eq. (PO]) . The second term in the exponent 
of the scale dependent factor on the right-hand-side of the above equation is caused by interactions in the last term 
of Eq. (|39p . Due to interactions each k mode is not frozen after horizon exit but evolves slowly. As larger modes 
spend less time outside the horizon they are less affected, which introduces additional k dependence. Both e and 
I77I in Eq. ([15|) are evaluated at horizon crossing. Although \r]\ ex. is a function of time this does not introduce 
additional scale dependence as all modes after horizon crossing are affected the same way by the evolution of \r]\. 
However, |7y| cx means that for cosmological scales |?7| is exponentially suppressed and the e term dominates 
the spectral tilt. In principle \ri\ introduces an anisotropic scale dependence. But as this term is subdominant, the 
direction dependence of \ri\ is suppressed. 

The presence of non-linear term in Eq. p9l) makes the vector field perturbation non-Gaussian. As we have chosen 
SWi to satisfy Gaussian statistics, non-Gaussianity is encapsulated in SW2- Thus SW2 satisfies the full non-linear 
equation. However, the curvature perturbation in the Universe is predominantly Gaussian. So if vector fields are to 
generate the dominant contribution to the curvature perturbation, they must be predominantly Gaussian too, that is 
< SWi- With this condition 6W2i can be seen as the second order perturbation. Thus perturbing Eq. ([39]) to 
the second order we find the equation of motion for (JVF^fj 



SW^i + 3HSW^, + 2V^, = 0. (46) 
Dropping out terms proportional to the slow-roll parameters in Eq. (j40p we write 

Assuming slow-roll holds we may also drop the first term in Eq. ()46p and write 

^^^^ = (48) 

Taking the Fourier transform it becomes 

5W^. (k) = - A I '^^0^5 (k - qi - q^) V^, (qi, q^) , (49) 

where V2i (qi, q2) is 

Vl\ (qi, q2, t) = jLi-bh^-dc [w^swl^ (q^) j^^rf (q^) + wfSWt^ {cii) 5W',, (q^) -|- WpW^, (qi) ^W^f, (qa)] ■ (50) 

As we are interested in the superhorizon evolution of the field perturbation, to find 5W2i we integrate Eq. P^)) from 
the horizon exit at tk^ where k/a itk) = 1, to some later time t. Because is slowly rolling and SW^^ is constant 
by definition with H k, const, the only time dependent term in Eq. ([501) is / oc a"**. Thus solving Eq. (|^^ we find 

5W^.(k,t) = -^ / ^^^<5(k-qi-q2)T^(k)y2';(qi,q2,0- (51) 

From this solution we can find the bound on the strength of self-coupling §1/ f for the condition 5W2 < SWi to be 

2 

consistent. Putting 6Wi ~ H into Eqs. ([50]) and (|5ip it follows < H. This ensures consistency of using second 
order perturbation theory to calculate 6W2 and that perturbations of vector fields are predominantly Gaussian. 



C. The Three- Point Correlation Function 



The three point correlation function from the classical evolution of the field is 

gf (ki, k2, kg) EE (^SW/ (ki) 5Wl (k2) SW!: (k3)) 



(52) 
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Because SWi ^ 5W2 the largest contribution to 173^ comes from the term of the form (^SWi5Wi^SW2) oc 
^ {6Wi5Wi5Wi -k 5Wi) , where a star denotes convolution. As one has to be careful in keeping track of indices 
we write the dominant term of 173^ explicitly 

at - I [{swl (ki) 6Wf^ (k2) SWt (kg)) + (swl (ki) SWl^ (k2) 6Wt (ks)) + {swl {k,)5Wf„, (k2) SWt (ka))" . 

(53) 

To evaluate this expression we use Eq. (jSip and Wick's theorem to express four point functions in terms of products 
of two point ones. After tedious algebra and using Eqs. we obtain 



gf = - (2.)^ S (k, + k. + k3) T^'itl''' (ki, k., ka) ^ ^^^^ 

where the anisotropy tensor 7^^)/^'' is defined in Eq. ([25]) and we dropped the gauge index from as all vector fields 
have the same spectrum. Taking V+ = (iJ/27r)^ we recover exactly the same result as obtained by a more tedious 
calculation in the quantum in-in formalism with the dominant term in Eq. (|29p . 



V. THE CURVATURE PERTURBATION 
A. SN Formula 

In the above we have calculated correlators of the field perturbation. However, field perturbation is not an observ- 
able, but the metric perturbation is. We choose a uniform density slicing in which the perturbation of the metric 
on superhorizon scales is described by the intrinsic curvature (. The easiest way to calculate ( is using the 6N 
formula |65, 6d|. This formula was first extended to include vector fields in Ref. and used for non-Abelian fields 
in Refs. [il Q 



C (x, t) = N^6<l, + N^SWt + -N^^ {5cj,Y + N^^M^ + ,^K^'5W^5W^ + . . . , (55) 



where 



t t ^ J 

N in these expressions is the number of e-foldings of local expansion from the initial flat hypersurface to the final 
uniform density hypersurface at final time t when ( becomes constant. In Refs. [21I, t was taken to be just after 
the horizon crossing. However, as Eq. ([29]) shows, the bispectrum of the field perturbation is actually dominated by 
the interaction of classical fields during classical evolution outside the horizon. 

Derivatives in Eq. (|56|) are taken with respect to homogeneous fields. The precise form of these derivatives depends 
on the mechanism through which the field perturbation generates (. S(f> in this equation is the perturbation of some 
scalar field if any of such fields contribute to the curvature perturbation. In this paper, by keeping only the second 
and fifth terms of the right-hand-side of Eq. (|55p we assume that predominantly vector fields contribute to ^ and any 
other source is negligible. 

B. The Spectrum 

Let us first consider the two point correlation function of In Fourier space we may write 

(C (k) C (k')) = NtN^ {SWt (k) 6W^ (k')) ^ {2nf S (k + k') |f T'c (k) , (57) 
where (5VFf = SW^, + \5W^^ + .... From Eq. ^ we find 

Vc (k) = ^ SabNtN'^Tf^ (k) V% [k) . (58) 

a.b 
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(k) i. 



Due to the presence of T,^ ( k ) in this expression the power spectrum of C, has an angular modulation. The isotropic 
part is 



(59) 



where the sum is over all gauge fields and Na is the absolute value Na = l-^f |- We also used the fact that power 
spectra 7^" = V+ are the same for all fields. Then the total spectrum of C is 



Vc (k) = V, 



ISO 

c 



Eg (N--k)' 



(60) 



If a = 1 this expression reduces to Eq. ([T]) with quadrupole anisotropy of an amplitude — —1 and such a large 
anisotropy is ruled out by observations. Thus, a massless U {!) vector field cannot generate the total C. However, 
with the large number of randomly oriented vector fields, the anisotropy is suppressed by the number of fields. To 
see this, note that N is proportional to vector fields VF". Since all vector fields satisfy the same equation of motion, 
assuming similar initial conditions it is reasonable to expect that contributions of all to N are of the same order 
In this case NaS are of the same order too, in particular Nw = Na ^ Nj, for all a and b. Then X^a-^a 
where Af is the number of fields and Eq. ([BU)) becomes 



Vc (k) = Afr+N'^ 



w 



1 - 



1 ^ 



w ■ k 



(61) 



where W° are unit vectors along the directions of homogeneous vector fields and we used the fact that W° = 
N° = 'N°'/Na- With Af vector fields, the anisotropic part of the spectrum is a sum of J\f quadrupoles. If these are 
randomly oriented, the anisotropy is suppressed by Af^^. Thus, with the large enough non-Abelian symmetry group 
the total curvature perturbation can be generated solely by vector fields. In view of contradicting conclusions of 
Refs. (39l . |40| we consider two bounds on 5^ to estimate Af. If we accept that a systematic effect causing the large 
anisotropy in the spectrum (0.29 ± 0.031 as claimed in Ref. [s^) is unknown, the bound on the anisotropy in of 
the primordial origin can be taken to be 17^ < 0.29. In such a case four vector fields A/" = 4 is enough for this bound 
to be satisfied. If, on the other hand, the large observed anisotropy in is caused by the WMAP beam asymmetry, 
as claimed in Ref. [i^, then the corrected bound on primordial statistical anisotropy gives jg^l < 0.07 [4lj. To satisfy 
this bound Af > 15 is needed. These estimates for Af are made assuming random orientation of the homogeneous 
vector fields and their similar magnitudes. If, however, all vector fields are parallel or one of the TV^'s is dominant, 
then = — 1 and such configuration is excluded. 

It is also possible to generate statistically isotropic curvature perturbation by considering a triad of orthogonal 
vector fields with equal norm. This configuration in a context of vector inflation was studied in Refs. |5l]. [52|. where 
SU (2) group is considered. However, we feel that a scenario with random orientation of larger number of fields is 
a more natural setup. In addition, such setup also has an advantage of providing observational signature, that is 
non-negligible statistical anisotropy. 



C. The Bispectrum 

In this section we calculate the bispectrum at the end of infiation. In momentum space it is defined by 

(C (ki) C (k2) C (ks)) = {27rf S (ki + k2 + ka) i^c (ki, ks, kg) . (62) 

Non-vanishing is the result of two contributions. The first contribution, let us denote it by _B^i. is from the 
non-Gaussian field perturbations due to self- interactions of the vector fields. The second, -6^2; is due to the non-linear 
terms in the SN formula in Eq. (|55p . 

Let us start by calculating the first contribution. From Eq. ([55]) 



(C (ki) C (k2) C (ka)) D n/n^N^: (^SW/ (ki) SWl (k,) SW^: (k,)) 



(63) 



13 



The three point correlation function of the field perturbation was calculated in sections IIIII and IIVI Taking the result 
in Eq. the bispectrum becomes 



Bci - -Send 



2 47r4 e 



^VlN^Nf^K^T^^^^^ (ki,k2,k3), 



(64) 



where ^ond = ffc/V/cnd is the strength of the self-coupling of canonically normalized vector fields at the end of inflation. 
Fallowing Refs. [17i. i45i] we introduce vectors A^" to simplify expressions for the bispectrum 



X^"(k)^7'+7VfI^Mk =V+Na 



k 



(65) 



where no summation over a is assumed. Using this definition and Eq. (|25p the bispectrum in Eq. (j64p becomes 



Ci 



-47r' 



fend 



c.p. 



(66) 



In this equation 'c.p.' stands for cyclic permutations of vectors k. 

The second contribution to the three-point correlator of the curvature perturbation is from non-linear terms in 
Eq. (j55p . As we assume that only vector fields generate Q the dominant term will be 



(C (ki) C (k2) C (ka)) D ^N^^N-^Nt {5W^ * 5W^ (ki) 6W^^ (ks) SW^ (kg)) + c.p. 
Using Wicks theorem and Eq. (1431) the bispectrum from the the above expression becomes 



(67) 



B^2 = 4^^ 



1 { 



f-M?(ki 



7V^.''X5(k3)-KC.p. 



(68) 



Both _B|^i and Bq2 depend not only on the absolute values of wavevectors k but also on their direction, making the 
bispectrum anisotropic. The total bispectrum is Bq ~ Bqi + Bq2- To evaluate which term is the dominant one note 
that Nab oc Na/W^, where Nat = l^jf]- Thus, in order for the first term to dominate, gcndW > H must be satisfied.^ 
If this is the case, the slow-roll condition \'rifj \ < 1 is violated. Then the evolution of the homogeneous modes of 
vector fields becomes strongly non-linear and the above calculations do not apply. However, if / is modulated by the 
inflaton or some degree of freedom which is stabilized at the end of inflation then f^nd is equal or very close to unity, 
i-e. /end ^ 1 and gcnd 9c- In this case we can expect gcW ~ iJ if is not much larger than H (see the footnote 
on page[S]), as it is natural for g,. to be not much bellow unity in particle physics models. If this is the case, then 
both contributions to Bq are comparable. If, on the other hand, gcndW < H then _B^2 contribution to the bispectrum 
dominates. 



VI. THE END-OF-INFLATION SCENARIO 



Let us implement the results of previous sections to a specific example using the end-of-inflation scenario. This 
scenario was suggested in Ref. [63 invoking only scalar fields. In usual hybrid inflation models inflation ends when 
the waterfall field is destabilized by the inflaton. This happens when the inflaton reaches some critical value and the 
waterfall fleld mass becomes tachyonic. As this critical value is determined solely by the inflaton itself, inflation ends 
on a uniform energy density slice. If, as suggested in Ref. [63|, this critical value is modulated by some additional fleld, 
then the uniform density slice does no longer coincide with the end-of-infiation slice. This induces the perturbation 
in the distance between fiat and uniform density slice, which is equal to the perturbation in In Ref. [l5| it was 
shown that if the modulating field is U (1) vector field, the generated C is in general statistically anisotropic (see also 
Ref. [17]). In this example we extend the scenario proposed in Ref. [l5i| to include non-Abelian vector fields. 



^ This is in contrast to Refs. |2ll . |22|| where the opposite bound was assumed. 
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A. The Model 



Let us consider a Lagrangian which is invariant under transformations of some non-Abehan symmetry group G 

L = \d^^d^v - \fFrF-^, + ixr [{D^<^>)^ D^'^] - V $) , (69) 

where Tr [. . .] stands for trace, is the inflaton field and is the field strength tensor defined in Eq. ([6]). The gauge 
kinetic function / may be the function of the inflaton / [ip). This has an advantage that we don't introduce additional 
degrees of freedom. The behavior of /((/?) for Abelian vector fields was studied in Refs. (32l. |33|. It was found that 
the required scaling / cx becomes an attractor solution in a large parameter space. Such kinetic function was 
also studied for non- Abelian vector fields in Ref. [ssj. However, for the present purpose we don't need to assume the 
source of modulation of /. 

$ in Eq. is the Higgs field corresponding to a non-trivial representation of G while the covariant derivative of 
the Higgs field is given by 

D^ = d^ + i\AT''Al, (70) 

where T° are generators which satisfy the Lie algebra [T°, T**] = if'^'"=T'^ of an unbroken symmetry group G and Aa 
is the gauge coupling constant, coupling the Higgs field to the vector gauge fields. 
The effective potential V in Eq. ([55)1 is taken to be 



V (</., $) = ^ [Tr + y^'Tr ($+$) +V(ip), (71) 

where A, k and M are constants with M being a symmetry breaking scale. V {(p) is the potential of ip providing the 
slow-roll inflation. V {ip, $) can be expressed in a more familiar form if we write the Higgs field as 

$ = 01, (72) 

where 1 is the matrix defining the direction of symmetry breaking in the field space with Tr [l^l] = 1. Then Eq. (j7ip 
becomes 



V {p, <f) = ^A {cf" - M^Y + i^V'^' + V{ip), (73) 

which is the potential of the hybrid inflation. But in contrast to the standard hybrid inflation scenarios we assume 
that the dominant part of the curvature perturbation is not generated by the inflaton field. Instead C, is generated by 
the gauge fields through the gauge coupling constants which couple them to the Higgs field in the covariant derivative 
in Eq. ([70|. To see this, note from Eqs. ([M)) . ((70)1 and ([75)1 that the effective mass squared of the Higgs field ^ is 



mis W = K V - AAf' - \\AIAI\^T:''tH. (74) 

In the unitary gauge 1 is such that the last term in this expression is diagonalised to obtain a sum of massive vector 
fields M"-^ A'^^A'^ . M"'' is a diagonal matrix with the only non-zero elements corresponding to broken generators. 

Note, however, that M°-^ is not the mass matrix of the vector fields but (fP'M"'^ is. Without the loss of generality, we 
can arrange generators T° in such a way that low a's correspond to generators of unbroken subgroup and higher a's 
correspond to the broken ones. Then M"** will have non-zero elements only in the lower right block, which can be 
written as M'^'^ , where we used bars over indices to remind us that they run only over the broken generators but not 
the full group. 



B. The Curvature Perturbation 



The curvature perturbation in this set up has two contributions 

C = + Cc- 



(75) 
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The first contribution C,^ is generated at the horizon crossing during the slow-roll inflation. The resulting power 
spectrum of is well known to be [68!] 



1 fHu^^ 



where iJ^. and are the Hubble and slow-roll parameters evaluated at the horizon exit and we used = dN/ dtp = 

(2mpjefc) for the slow-roll inflation. The contribution of to non-Gaussianity is proportional to slow-roll pa- 
rameters at horizon exit and, therefore, too small to ever be observable (ssl. [69|. 

When the inflaton crosses some critical value Lp^, at which rn^g in Eq. (|74p becomes negative, inflation terminates 
and the Higgs field rolls down to the minimum of the potential. From Eq. (|74p we find 



pI = 4^^' - -^W^WfW^ (77) 

K K J 



where we made use of the temporal gauge and Eq. ([5]) to specify tp^ in terms of the physical, canonically normalized 
vector fields W"". The second term in Eq. ([77)) is subdominant, i.e. k^pI ~ \A-P , giving 

AAf2 > M'^^Wfwf/f. (78) 

But, due to perturbations of the vector fields, it modulates the critical value of the infiaton, making p^ (x) a function 
of space coordinates x. Thus the end of inflation hypersurface does not coincide with the uniform energy density 
hypersurface which results in the generation of the curvature perturbation C,c- Up to the second order C,^ is given by 



Cc = N,6p,+N,,{5p,f , (79) 



where iVo = dN/dipc = {2mp^ee) , Ncc = d'^N/dipl and €c is the first slow-roll parameter at the end of infiation. 
The perturbation of pc can be written as 

where Sf = f Sp^/'^c is the variation of the kinetic function / corresponding to the time shift from the hypersurface 
of the uniform energy density to the end of inflation. As we require / c>c a"^ for the gauge and self-couplings of 
the canonically normalized fields to be small, the time derivative of / is negative, / < 0. Thus, the second term in 
Eq. (j80p suppresses Cc- Even more so, if this term dominates Spc, no perturbation is generated at the end of infiation. 
To ensure, this does not happen, we require the first term to dominate, which gives the constraint 



K'^feipcmpi I 



« ec, (81) 



where /, is evaluated just before the end of infiation. To simplify calculations we assumed a stronger condition, that 
the second term in Eq. (j80p is completely negligible. Then Sipc is equal to 

Jefc 

and the isotropic part of the power spectrum of Cc in Eq. (|59[) is given by 

nr=^+E^', (83) 



where from Eqs. (|79p and 



M°-''W'' 

N^ = -N-^r^- (84) 



^ Note that the second term in this equation was neglected in Refs. Il5l llTll . 
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As was mentioned in subsection |VB] if only one vector field contributes to the curvature perturbation, the anisotropy 
in the power spectrum is — and such a large value is excluded by observations. This is the case, for example, 
with the Abelian vector field. Thus authors of Refs. (isl. ITtI assumed that Co < which case V^^° is dominated 

by the scalar field contribution and the subdominant vector field contribution generates anisotropy in the spectrum 
of C with \g(^\ < 1. In our case, since we are dealing with non-Abelian vector fields, not one but several vector fields 
contribute to C- If their orientation in space is random, the anisotropy in the spectrum is suppressed by the number 
of fields J\f (see Eq. (|5T|) ). Thus with the large enough J\f (which is evaluated in subsection I VI C[) we can generate 
the total curvature perturbation without violating observational bounds on g^. V(^^ dominates the spectrum of the 
curvature perturbation if jA^al ^ N,^. Using Eq. ([Mil this bound becomes 



.^^)'»- = e-'^"''' (85) 



W 

where X\/ fc is the gauge coupling of the vector field to the Higgs field at the end of infiation. Evaluating Eq. 
we assumed that all gauge fields are of the same order, i.e. W ^ W"" for all a. We also used the fact that absolute 

values of matrix elements of generators T° in Eq. ([7(I| are of order unity so that Tr (^M"^^) - A^. The slow-roll 

parameter 77 in Eq. (1851) is 77 = mpjV^^,^ ((/j) /V {(p) and subscripts denote the second derivative of V {(p) with respect (p. 
As discussed in Ref. [63| with Ce dominant, 77 has nothing to do with the spectral index of the curvature perturbation 
and can even be 77 ~ 1. In this case the right hand side of Eq. ([55)) can be far bellow unity. 

Comparing two bounds in Eq. ([5T|) and (|85l) we find that for successful end-of-inflation scenario, in which gauge 
fields generate the dominant contribution to the curvature perturbation, the homogeneous gauge field value must 
satisfy 



From the first and last terms we find 



\mpij \mpij 

where the first constraint is explained in the footnote on page HI Assuming the scale of inflation to be of the order 
of GUT scale, i.e. ^ 10^^ GeV, H/mp\ ~ 10~^. Taking ek, when cosmological scales leave the horizon, to be of order 
10-2 , the bound in Eq. ^ gives IQ-* < W/mpi < 10"^ 



C. Anisotropic Spectrum and Bispectrum 



To find the full power spectrum of the curvature perturbation generated by the non-Abelian gauge fields let us 
substitute Eq. ^ into Using Eq. we find 



Vc (k) = V+C^ {M^f" W^'W 



(M2)'''(w"-k) (W^.k) 



where (M^)" ee M^^M^^ is a diagonal matrix, W is the modulus of the vector field W"' = \ \ and C is defined as 

1 



C . 



(89) 



Note that Vq^ is determined solely by the massive vector fields. If the homogeneous values of all vector fields are of 
the same order, i.e. W ~ W"' for all a, then the power spectrum in Eq. ([88)) becomes 



(k) « \\MV+ {CWf 



(90) 



where M is the number of massive vector fields and we used Tr was discussed after Eq. ([6T)) A/" > 15 

or A/" > 4 is needed to avoid observational constraints on , depending if the systematics causing detected anisotropy 
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in the spectrum is believed to be the asymmetry of WMAP beams or unknown. The Lagrangian in Eq. (j69p was 
assumed to be invariant under the transformation of some non-Abehan symmetry group G. To estimate the minimal 
rank of the group G, which satisfies the bounds on A/", let us assume that G is a special unitary group SU [N) which 
at the phase transition is broken to SU {N — 1). Such symmetry breaking results in A/" = 2N — 1 massive gauge fields. 
Thus for a weaker bound on JV the SU (3) group already generates < 0.29. If, on the other hand, the stronger 
bound on Af applies (with < 0.07), at least SU (8) is needed. In realistic particle physics models one has to be 
careful in choosing a symmetry group as not to overproduce monopoles after the symmetry breaking f70'|. However, 
the above estimate demonstrates that the primordial curvature perturbation can be generated with gauge fields of 
reasonably large groups. Even more so, the anisotropy in the spectrum g^, generated by such groups, are of the 
magnitude which will be possible to test in the very near future by the Planck satellite. As it is shown in Refs. (36l.l4T| 
the Planck data will allow to constrain with an accuracy of 0.01. 

To find the bispectrum for the non-Abelian end-of-inflation scenario let us differentiate N°- in Eq. one more 
time 

Nt]' = -S,,CM^'' + ^-^. (91) 

Using Eq. (j78p and K^ifi'^ ~ XM^ one can easily check that the first term dominates in this expression. 

Also let us find the vector Mf (k) introduced in Eq. ([55]). With Nf- calculated in Eq. Ml (k) becomes 



(k) = -V+CAr^ Wf - k, (w^ ■ k) ■ (92) 



The first bispectrum Bqi in Eq. which is due to self interactions of gauge fields, can be calculated using the 

results in Eqs. (|84)) and (|52)) . After some algebra we obtain 



Bi^i = Att^^^^^^C^VI (f''''''f''^^ + r3hfabf\ ms'^M^'^M''^ {W'' ■ W 



W^' • - 2 (W^ ■ ki) ■ ki) + fki • ks) ■ kil f W • kg) + c.p.j . (93) 

First, note that the bispectrum is anisotropic due to its dependence on the direction of k wavevectors. The anisotropy 
in the bispectrum from self-interactions is solely determined by massive vector fields. The amplitude of _B^i however 
is determined by vector fields and structure constants of the whole group as it has unbarred indices. 

To evaluate the magnitude of B^i let us assume that all vector fields are of the same order W and Tr ^Af''^ ~ X\. 

Using the isotropic part of the spectrum in Eq. (|90p and assuming that structure constants are of order unity we find 

~ > I ) rpf- ^^^^ 

The expression for the second part of the bispectrum B(^2 is quite long and the full result is given in Appendix |B] 
But it is easy to see from Eqs. (jHH), ([Ml) and ([55]) that it is anisotropic too and it is determined solely by massive 
vector fields. The isotropic part of Bq2 is 

Note, that compared to the single field end-of-inflation scenario in Ref. [b^I, the bispectrum is suppressed by the 
number of fields. However, as discussed after Eq. (|90p we do not expect Af to be too large. Also note, that the 
bispectrum from the self-interactions B(^i has an additional factor {gcndW/ H) , where gcnd is the strength of self- 
coupling of the canonically normalized gauge fields at the end of inflation. Although this factor can not be much 
larger than one, as it would make the evolution of gauge fields strongly non-linear, it might be not much smaller 
than unity. Finally as is shown in Eq. (jB2[) the anisotropy in the bispectrum as in the spectrum is suppressed by the 
number of massive gauge vector fields. 



VII. SUMMARY AND CONCLUSIONS 



The possibility of vector fields playing a non-negligible role in the very early Universe is attracting more and 
more attention both from theorists as well as data analists. The role of vector fields can be to provide either an 
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anisotropically expanding Universe during inflation or directly affecting or even generating the primordial curvature 
perturbation or both. Some works study the possibility of vector fields driving inflation too. Because vector field, in 
contrast to a scalar field, chooses a preferred direction, the smoking gun of such models is the statistically anisotropic 
curvature perturbation. The effects of such anisotropy can be observed in temperature and polarization irregularities 
of the CMB sky. Indeed, with the measurements of Planck satellite, which is currently collecting data, it will be 
possible to constrain statistical anisotropy at the level of 0.01 [4l|. 

In this paper we studied the curvature perturbation generated by non-Abelian vector fields. Non-Abelian vector 
fields are one of the main building blocks of the standard model of particle physics and indeed of any gauge theory 
and their existence is an experimentally confirmed fact. Moreover, theories beyond the standard model contain 
large numbers of such fields. In this paper we consider massless non-Abelian fields with the Lagrangian of the form 
£ ^ -jf it)F^''F^^, with F^^ = df.A'^ ~ d^A^ + gj"'"' A'^^Af, and i"'"' being structure constants of a general Lie 
group. (7c is normalized in such a way that / (ig) = 1 when it is stabilized at some time tg. The perturbation spectrum 
for Abelian vector fields is fiat if the kinetic function scale s as / oc a~^^^. In the non-Abelian case however, / cx 
corresponds to a strong effective self-coupling g (t) = gd \/ f (t) for the physical vector fields — ^/jAf/a, which 
results in strongly non-linear evolution. For this reason we assume / oc a~^, which is also necessary for the field 
perturbation to become classical. The requirement for predominantly linear evolution of and its perturbation also 
puts the bound on the self-coupling of fields. As strong non-linearity affects all scales (not only the ones leaving the 
horizon) the bound 5c < 1 must be satisfied. 

In Refs. [32, 3^ it was shown that the scaling of the form / oc a can be achieved dynamically through the 
backreaction of vector fields on the evolution of the infiaton, if / is modulated by the inflation and the vector fields 
are Abelian. However, this induces anisotropic expansion of order the slow-roll parameter e. Anisotropic expansion 
in its own right introduces an additional source of statistical anisotropy in the curvature perturbation [26]. However, 
in our analysis, we do not require / to be necessarily modulated by the infiaton. To avoid anisotropic expansion we 
also assume a negligible contribution of the vector fields to the the overall energy budget during inflation. 

With this setup in section iHll we calculate the bispectrum of the fleld perturbation resulting from interactions of 
flelds. To calculate the three point correlation function at the tree level we employ the full quantum perturbation 
formalism, the so called "in-in formalism". The interaction Hamiltonian with the above Lagrangian has two terms, 
the cubic term with derivative couplings and the quartic term. The contribution from the first one is suppressed by 
the physical momentum p <C as compared to the second term. While the bispectrum from the quartic term is 
dominated by the classical evolution of fields, i.e. by interactions after a mode exits the horizon. 

This being the case, it is much easier to calculate the correlation functions from the homogeneous classical equation 
of motion. Such calculation is performed in section HVl It is shown that the result from this method is indeed exactly 
equal to the dominant part from the full calculation using the in-in formalism. 

In section IVl we calculate the spectrum and bispectrum of the curvature perturbation. It is found that the spectrum 
has angular modulation. However, in contrast to the single vector field case, the anisotropy g^; in the spectrum is 
suppressed by the number of fields (assuming random orientation). Thus, reasonably large groups can generate small 
but observable g^^. In Ref. fZTl it is shown that with the Planck data it will be possible to constrain g^ with the 
precision up to 0.01. 

The bispectrum of the curvature perturbation has two contributions, one from the non-Gaussian field perturbation, 
_B^i, and the other from non-linearity in generating -6^2- Both of those contributions have an angular modulation 
and are comparable if gcndW/H ^ 1, where W is the modulus of the homogeneous part of the vector fields and gcnd 
is the self-coupling strength of canonically normalized vector fields at the end of infiation. If this ratio, however, is 
much larger than 1, the evolution of vector fields is strongly non-linear and above calculations do not apply. In the 
opposite regime _B^2 dominates. 

In the last section I VII we present an example of a mechanism for vector fields to generate C. In this example we 
consider a scenario in which the curvature perturbation is generated by varying gauge coupling(s), which couple the 
Higgs field to vector bosons in the covariant derivative. In such models ( is generated by the vector boson fields 
corresponding to broken generators after the phase transition. We calculate the spectrum and the bispectrum. The 
anisotropy in both of them is suppressed by the number of massive vector gauge bosons. 
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Appendix A: Calculation of Integrals Z'*^ and I^^^ 



111 this appendix we show how to calculate integrals in Eqs. p6| and ((33|) . Let us rewrite here the first integral 



= ^Re 
Jo 



dr' (1 - ikiT') (1 - ik2T') (1 - iksr') e**^'^' 



(Al) 



During quasi-de Sitter inflation with H sa constant the scale factor is a w —1/tH. Using this and denoting x = k^r, 
the above integral can be rewritten as 



Jo 



smx cos a; smx cosa: , , 

Al A2 — I dx, 



(A2) 



where xq = ktTo — — oo corresponds to the initial time when modes are deep within the horizon and Xend = fctTcnd 
is at the end of inflation. Note that assuming all three fc's leave the horizon at similar time, Nk — — ln(|a;ond|) is the 
number of e- folds from the horizon crossing to the end of inflation. Ki and K2 are defined in Eq. ((28)) . 

The total integral in Eq. (jA2p is the superposition of integrals J sinx/x" dx and J cosx/x" dx, with n being a natural 
number. The order of n within each integral can be reduced integrating by parts until we arrive at superposition of 
terms sin Xcnd /x",j(j and cos Xcnd/x^j^^ with appropriate constants and the integral J cosx/xdx. The last one can be 
evaluated as follows. Let us write 



cosx 



dx = 



1 — cos X 



dx 



cosx 



dx 



din: 



(A3) 



Taking the limit Xend — >■ —00 and xq — first two terms in Eq. (|A3[) are equal to the Euler-Mascheroni's constant 
7 ~ 0.577 and the last term is —Nk- Expanding the result around Xcnd — ^ and neglecting terms proportional to 
x",jjj with n > 0, we arrive at the final expression in Eq. ([?7)) . 

The same method can be used to evaluate I^-^K The difference is, that integrating by parts the lowest order integral 
becomes 



sm X TT 

dx = -, 

X 2 



(A4) 



where we have taken a limit xq 



-00 and Xcnd 0. 



Appendix B: The Bispectrum B, 



C2 



The expression for the second part of the bispectrum B^2 can be calculated using the first, dominant term of N°j' 



in Eq. ([91 



Mt iki)NtJ'M'; (ks) = -VlC^ {M^f^ [(W^^ • W^") - (W^^ • ki) (w^" • ki) - 

• kg) (W^^ • kg) + (ki • kg) (1^' • ki) [W" ■ kg) 



where (M'^)'"* = M^^'M^^Af^'' is the diagonal matrix. Inserting this result into Eq. gives 



(Bl) 



C2 



47r^ 



■^ISO 



kl (w°-ki)^-/cf (w^-kg) 



^k\ (kl • kg) (W^ • kl) (W"^ • kg) + C.p.] lY^kiy 



(B2) 



where we assumed W ~ W"" for all a and Tr 



(A/3) 



ah 



A^. Note that the anisotropy in the bispectrum is suppressed 



by the number of fields, the same suppression as in the spectrum. 
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